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Abstract. The relation between the spectral decomposition of a self-adjoint operator 
which is realizable as a higher order recurrence operator and matrix-valued orthogonal 
polynomials is investigated. A general construction of such operators from scalar-valued 
orthogonal polynomials is presented. Two examples of matrix-valued orthogonal poly- 
nomials with explicit orthogonality relations and three-term recurrence relation are pre- 
sented, which both can be considered as 2 x 2-matrix-valued analogues of subfamilies of 
Askey- Wilson polynomials. 



1. Introduction 

Matrix-valued orthogonal polynomials date back to the 1950's in the work of M.G. Krein, 
see e.g. references in [2], [3]. More recently, matrix- valued orthogonal polynomials are stud- 
ied from an analytic point of view. In particular, analogues of many classical results in 
the theory of ordinary (scalar-valued) orthogonal polynomials have been generalized to 
the situation of the matrix-valued orthogonal polynomials, such as e.g. the three-term 
recurrence relation, the spectral theorem (Favard), theorems of Markov, Blumenthal, etc., 
see the overviews [2], [3] and references given there. Many examples of the general theory 
of matrix- valued orthogonal polynomials are motivated by matrix- valued differential equa- 
tions, see also [8]. Some of these examples are motivated from the well-known families of 
orthogonal polynomials in the Askey scheme |l4j . so the matrix- valued weight function is 
given by the scalar weight function times a suitable matrix-valued function. So in this case 
matrix-valued analogues of classical orthogonal polynomials, such as Jacobi, Laguerre and 
Hermite polynomials, are obtained. This theory sofar gives matrix- valued analogues of hy- 
pergeometric orthogonal polynomials. Very little is known about matrix-valued analogues 
of g-orthogonal polynomials. 

Another way of obtaining matrix-valued orthogonal polynomials is from group theory 
using matrix-valued spherical functions. An important case study has been given by 
Grunbaum, Pacharoni and Tirao [7j, in which they obtain matrix- valued orthogonal poly- 
nomials from the symmetric pair (SU(3),U(2)) by studying eigenfunctions to invariant 
matrix-valued differential operators. Again these matrix-valued orthogonal polynomials 
are analogues of a subfamily of Jacobi polynomials. In [15], [16] a different approach 
to such a group-theoretic approach has led to matrix-valued Chebyshev polynomials in- 
cluding relevant group theoretic interpretations of the construction, the three-term re- 
currence relation, weight function, differential equations, etc, using the symmetric pair 
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(SU(2) x SU(2), SU(2)). Again, in these cases the weight function resembles the cor- 
responding scalar weight function times a suitable matrix-valued function. Again, no 
g-matrix-valued orthogonal polynomials have yet emerged from this approach. 

In this paper we discuss a new way to obtain matrix-valued orthogonal polynomials with 
an explicit three-term recurrence relation as well as explicit orthogonality relations. In the 
examples it is clear that the weight function is not of the form of a classical weight function 
times a matrix-valued function. The idea is to look for the spectral decomposition of a 
self-adjoint operator which can also be realized as a higher order recurrence operator. In 
order to motivate the construction, we first note that if we consider an operator which can 
be realized as a 2N + 1-recurrence operator, the case N = corresponds to eigenfunctions. 
The case N = 1 is the case of the J-matrix (or tridiagonalization) method, which is used in 
physics to determine the spectrum of certain physically relevant operators, see [TO] , [12] and 
references given there. In [11] a more general method to obtain suitable tridiagonalizable 
operators is discussed. In this paper we restrict ourselved to self-adjoint operators that 
can be realized as 5-term recurrence operators and for which we have an explicit spectral 
decomposition. We show in Theorem 12.11 how this gives rise to 2 x 2-matrix-valued or- 
thogonal polynomials with an explicit (matrix-valued) three-term recurrence relation and 
explicit matrix-valued orthogonality relations. Because of computability reasons we stick 
to the 2 x 2-case, but we expect that it is possible to extend to larger size matrices. In 
Section H] we discuss an explicit example with an easy matrix- valued three-term recurrence 
relation, but an involved, but explicit, expression for the matrix- valued weight function. 
In Section [3] we discuss a general set-up, which is motivated by [TT] . and we work out a 
specific example in Section 1X21 which is related to the example in [TTj §4]. This motivates 
us to view the family of matrix-valued orthogonal polynomials discussed in the example of 
Section 13.21 as analogues of a subfamily of Askey- Wilson polynomials. 

As is well-known, it is very hard in general to obtain explicit expressions for the orthog- 
onality measures or weights for orthogonal polynomials defined by a three-term recurrence 
relation. The cases of associated classical orthogonal polynomials (in the Askey-scheme 
[T4"] ) amply demonstrates this point, see e.g. [13J for the case of two families of the as- 
sociated Askey- Wilson polynomials. It is therefore remarkable that we can obtain in this 
setting an explicit, even though complicated, expression for both the weight function and 
the three-term recurrence relations for the 2 x 2-matrix-valued orthogonal polynomials in 
the examples considered in this paper. Moreover, to our best knowledge this is the first 
instance of matrix-valued orthogonal polynomials that can be considered as matrix- valued 
orthogonal polynomials in a yet-unknown (possible) g-scheme of matrix-valued orthogonal 
polynomials, see [H] for the scalar case. Note that we do not have explicit expressions for 
the 2 x 2-matrix-valued orthogonal polynomials, and it would be of interest to obtain such 
expressions for these polynomials in terms of (yet to be developed) matrix-valued basic 
hypergeometric series of higher type, see Tirao [19] for the matrix-valued analogue of the 
classical hypergeometric function. 
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2. Matrix-valued orthogonal polynomials from 5-term operators 

In this section we study the relation between a self-adjoint operator realizable as 5-term 
operator and corresponding 2 x 2-matrix-valued orthogonal polynomials. The three-term 
matrix- valued recurrence relations for these polynomials follow from this realization of the 
operator, whereas the orthogonality relations for these polynomials follow from the spectral 
decomposition of the operator. The precise relation is given in Theorem 12.11 

We assume that we have an operator T on a Hilbert space H of functions. For T we 
typically consider a second-order difference or differential operator. We assume that T has 
the following properties; 

(a) T is (a possibly unbounded) self-adjoint operator on H (with domain D in case T is 
unbounded); 

(b) there exists an orthonormal basis {/ n K£Lo of H so that f n G D in case T is unbounded 
and so that there exist sequences (a n )^ =0 , (& n )J£L > ( c n)n°=o °f numbers with a n > 0, 
c n G R, for all n G N so that 

T f n = O'nfn+2 + &n/n+l + C n f n + 6 n _i/ n _i + a n - 2 fn-2- (2-1) 

In (b) we follow the convention that a_i = a_2 = 6-1 = 0. 

Next we assume that we have a suitable spectral decomposition of T. We assume that 
the spectrum is simple or at most of multiplicity 2, and we leave it to the reader to extend to 
higher order spectra. We assume that the double spectrum is contained in ^2 C cr(T) C R, 
and the simple spectrum is contained in Qi = a{T) \fl2CK. Consider functions / defined 
on er(T) CKso that /|q x : Qi — > C and /|n 2 : Q 2 C 2 . We let a be a Borel measure on 
Qi and V p a 2 x 2-matrix-valued measure on f2 2 as in [2J §1.2], so F: ^2 ~^ M 2 (C) maps 
into the positive semi-definite matrices and p is a positive Borel measure on Q 2 - Next we 
consider the weighted Hilbert space L 2 (V) of such functions for which 

/ \f(X)\ 2 da(X) + [ f*(\)V(\)f(\)dp(\) < 00 

and we obtain L 2 (V) by modding out by the functions of norm zero. The inner product is 
given by 

(f,9)= ! f{X)W)da{\) + [ g*(\)V(\)f(\)dp(\). 
J fix Jn 2 

The final assumption is then 

(c) there exists a unitary map U : % —> L 2 (V) so that UT = MU, where M is the 
multiplication operator on L 2 (V). 

Under the assumptions (a), (b), (c) we link the spectral measure to an orthogonality 
measure for matrix- valued orthogonal polynomials. Apply U to the 5-term expression (12. ip 
for T on the basis {/„}^L , so that 

\{Uf n ){\) = a n {Uf n+2 ){\) + b n (Uf n+1 )(X) 

+ c n (Uf n )(X) + K^(Uf n -t)(X) + a n _ 2 ([// n _ 2 )(A) (2.2) 
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to be interpreted as an identity in L 2 (V). Restricted to Qi (12. 2p is a scalar identity, and 
restricted to Q 2 the components of Uf(X) = (Uif(X), U 2 f(X)Y satisfy fFO) . 

In general, a 2N + 1-term recurrence relation can be solved using N x iV-matrix-valued 
orthogonal polynomials, see Duran and Van Assche [I]. Working out the details for N = 2, 
we see that we have to generate the 2 x 2-matrix-valued polynomials by 

\P n {\) = A n P n+1 {\) + B n P n {\) + A* n _ x P n ^{X), 

. a 2n \ = (c* b 2n \ (2.3) 



b 2n jf.\ a 2n+ \ J ' \p2n C 2n+ i 

with initial conditions P_i(A) = and Pq(X) is a constant non-singular matrix, which we 
take to be the identity, so Pq(X) = I. Note that A n is a non-singular matrix and B n is a 
Hermitian matrix for all n G N. Then the C 2 -valued functions 

Un{X) = ( T Yf*M) , Z&A) = (r^/^l) , Ul{X) - ( U *f»M 



satisfy (12.31) for vectors for A G Ox in the first case and for A G Q 2 in the last cases. Hence, 

Z4(A) = P„(A)W (A), Z^(A) = Pn(A)W 1 (A), W 2 (A) = P n (A)W 2 (A), (2.4) 

where the first holds cr-a.e. and the last two hold p-a.e. We can now state the orthogonality 
relations for the matrix-valued orthogonal polynomials. 

Theorem 2.1. With the assumptions (a), (b), (c) as given above, the 2 x 2-matrix-valued 
polynomials P n generated by (12.31) and P_ 1 (X) = 0, Po(X) = I satisfy 

P n (X)W 1 (X)P m (X)*da(X) + [ P n {X)W 2 {X)P m {Xydp{X)=8 nm I 

JClo 



where 



and 



W (X) - ( |/7/o (A)l Uf (X)Ufi(X) ] 

WliX) - [ umu fl (x) \u fl (x)\> »' °~ a - e - 



W (\) f( U h( X )> Ufo(X)) VW (Ufo(X), Uh(X)) v(x) \ n _ np 
W2[A) ~ \(Uf 1 (X) : Uf (X)) v{x) (?7/ 1 (A) > l7/i(A)>v W> /' p ^ 

where (x,y) v{X ) =x*V{X)y. 

Theorem 12. II can be phrased more compactly, and then the generalization to self-adjoint 
operators T realizable as higher order recurrence relations can be phrased compactly as 
well. Since we stick to the situation with the assumptions (a), (b), (c), the multiplicity of T 
cannot be higher than 2. Note that the matrices W\(X) and W 2 (X) are Gram matrices. In 
particular, det(Wi(A)) = for all A. So the weight matrix Wi(A) is semi-definite positive 
with eigenvalues and tr(W(A)) = \Uf Q {X)\ 2 + \U h(X)\ 2 > 0. Note that 



ker(W(X)) = C( U S^) = Si!! , ker(H/(A)-tr(iy(A))) = C' f//o(A) 



-uf (x)J \ufi(x)J ' "v"vw -yuMX) 

Moreover, det(W2(A)) = if and only if Uf (X) and Ufi(X) are multiples of each other. 
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^ /I 0\ / (/ 2n , f2m)n \f2n1 hm+l/n 
nm \0 lj \(f 2n +l, f2m)n (fln+U flm+l/Hy ^ r> _ , 

/ (U f 2n ,U f 2m ) L 2(y) (U f2n,U f 2 m+l) L 2 (V) 
\(^ /2n+l, ^ /2m)L 2 (V) /2n+l, ^ }2rn+l)v 2 {V) / 



Split each of the inner products on the right hand side of (12. 5p as a sum over two integrals, 
one over f2i and the other over Q 2 . First the integral over Vt\ equals 



4 Uf 2n {X)UUJ^)da{\) J Qi Uf 2n (X) Uf 2m+l (X)d a(X) \ 
J ni Uf 2n+1 (X)Uf 2m (\) dcr(X) J Qi Uf 2n+1 (X)Uf 2m+1 (X) da(X) J 



Uf 2n {X)Uf 2m {X) Uf 2n {X)Uf 2m+l {X) 



da(X) 



Ui \Uf 2n+ i(X)Uf 2m (X) Uf 2n+1 (X)Uf 2m+1 (X) / 
P n (A)W 1 (A)P w (A)*da(A), 



where we have used (12.41) . For the integral over Q 2 we write Uf(X) = (Uif(X),U 2 f(X)Y 
and V(A) = (%(A))f,- =1 , so that the integral over fi 2 can be written as 



2 



0,2 



V / [ U jf2nW v ijW U ihmW Ujf2n(X)v ij (X) U i f 2m+1 (X) \ , . 

^ ^,/2n+i(A)^(A)f/ 4 / 2m (A) ^•/ 2n+1 (A)^(A)^/ 2m+1 (A)y' ^ l J 



V /" P(x\( u M x A ( u M x )Xp 

— 1 



(A)%-(A) dp(A) 



Pn(A)iy 2 (A)P m (A)*rfp(A), 



'n 2 
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where we have used (12. 4p again and with 

w(x) - V ( U M*)\ ( u M x Kv (X) 



2 



V,; m (UjM WfoW UjM WfiM ) (2.8) 
~ ^ ] WiW^/o(A) U j fx(X)U i f 1 (X)) 



f(Uf (X))*V(X)Uf (X) («7/i(A))*y(A)l7/ (A) 
■ ^(C//o(A))*y(A)i7/ 1 (A) (17/i(A))*V(A)17/i(A) 

and putting (j22j) and ([22]), ([23]) into ([23]) proves the result. □ 

In case we additionally assume T is bounded, so that the measures a and p have compact 
support, the coefficients in (12. ip and (I2.3P are bounded. In this case the corresponding 
moment problem is determinate, see [2] Thm. 2.11], and Theorem 12.11 gives the explicit 
expression for the weight function. 

Remark 2.2. Assume that Oi = cr(T) or ^2 = 0, so that T has simple spectrum. Then 

C 2 (W 1 da) = {/: R^C 2 | / f(X)*W 1 (X)f(X)da(X) < oo} (2.9) 
has the subspace of null-vectors 



Af = {/ G £ 2 (M^<x) | / 



/(A)*W(A)/(A) cfo-(A) = 0} 



{/ G £ 2 (WW | /(A) = c(A) ( ) a _ a . e .}, 



where c is a scalar-valued function. In this case L 2 (V) = C 2 {W\da) / N . Note that U n : R — » 
L 2 {W\da) is completely determined by Ufo(X), which is a restatement of T having simple 
spectrum. From Theorem 12.11 we see that 

{P n {-)vx, P m {-)v2) L 2 (Wida) = 8 nm {vi,v 2 ) 
so that {-Pn(-) e i}ie{i,2},neN is linearly independent in L 2 (Wider) for any basis {ei, of C 2 . 

3. A GENERAL CLASS OF EXAMPLES 

In [TTJ we have studied a general procedure to obtain self-adjoint tridiagonalizable op- 
erators, and in this section we show how to extend this to obtain self-adjoint operators 
which can be realized as 5-term recurrence. This brings us back to the situation of Section 
12] hence leading to 2 x 2-matrix- valued orthogonal polynomials. Of course, we still need to 
obtain the spectral decomposition of such operators as well. We extend [TTJ §2] in Section 
13.11 and we present an example of the construction using little g-Jacobi polynomials in 
Section 13.21 The analogue of the Jacobi polynomials is rather involved, in particular the 
spectral decomposition, and this is worked out in [5J. 
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3.1. Self-adjoint penta-diagonalisable operators. Let p, and v be positive Borel mea- 
sures with finite moments on the real line R so that p is absolutely continuous with respect 
to v. Let r = ^ be the Radon-Nikodym derivative, so r > 0. We assume that we have 
a (possibly unbounded) self-adjoint operator L on L 2 (p) preserving the space of polyno- 
mials in L 2 (p) and the existence of an ortho normal basis {3> n }neN of L 2 (p) of polynomial 
eigenf unctions of L, so L$ n = A n $ n , A n G R. Moreover, we assume the existence of an 
orthonormal basis {<fi n } c £ =Q of polynomials in of L 2 [y) such that for all n G N 

(f) n OL n $ n -|- f3 n $ n _i /3 n , 7n GR (3.1) 

(with the convention /3 = 7o = 7i = 0). We assume that the polynomials are dense in 
L 2 (p) and L 2 [y). Finally we assume that the Radon-Nikodym derivative r is a polynomial, 
necessarily at most of degree 2 by (13.11) . We denote by M(r) and M(x) the multiplication 
operator by r and by x. 

Lemma 3.1. T p = M(r)(L + p) , p G R, is a symmetric five-diagonal operator on L 2 {y) 
with respect to the orthonormal basis {0 n }^L o ; 

T p 4> n = a n (j) n+ 2 + b n (j) n+1 + c n (f) n + b n -x(j) n -i + a„_ 2 0„-2 

where 

a n = a„7n+2(A„ + p), b n = a n /3 n+ i(X n + p) + /3 n (A n -i + p)ln+u 
c n = a 2 n (\ n + p) + /3^(A n _i + p) + 7n(-^n-2 + p)- 

Proof. This is completely analogous to [Tij §2.1]. □ 

Since the orthonormal basis {4> n }^=o of L 2 (v) consists of polynomials, we have 

x(j) n (x) = 8 n (f) n+1 (x) + £„0„(x) + 6»„_i0„_i(x), (3.2) 

for 9 n , G 1, S n ^ for all n G N and the convention #_i = 0. 

Corollary 3.2. T^ r = M(r)(L + p) + tM(x), p, t G R, is a symmetric five-diagonal 
operator on L 2 (v) with respect to the orthonormal basis {0 n }^=o? 

T P ' T (f) n = a n (f) n+ 2 + b n (f) n+1 + C n (j) n + bn^n-x + a n _ 2 0„-2 

where b n = b n + r^ n , c n = c n + r^ n; and £/ie notation as in Lemma \3. 11 

Note that in case L is a second-order differential or difference operator, then so is T. 
However, the coefficients of T get more complicated and in order to carry through the 
programme of Section [2] we need to be able to calculate the spectral decomposition of T P ' T 
for suitable p, t as well in another way. 

Remark 3.3. It is clear that we can extend this to higher order recurrences. So if we 
assume r to be a polynomial of degree N and the recursion (13.11) to have iV + 1 terms, we 
end with a 2N + 1-recursion for the operator in Lemma 13.11 and Corollary 13.21 
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3.2. Example: case of little g-Jacobi polynomials. We work out the details of the 
general programme of Section 13.11 for the case of the little q- Jacobi polynomials, cf. [Hi 
§4]. Let, as usual, < q < 1, and we follow standard notation for basic hypergeometric 
series as in [5], see also (9], [13] . 
The little q- Jacobi polynomials are 

p n (x) = p n (x;a,b;q) = 2 <Pi ! Q, qx ) (3.3) 

V aq J 

with leading coefficient 

l n (a,b) = (-1)^-1^-1) (3 . 4) 

(a?; Q)n 

and for < a < q^ 1 , b < q^ 1 the little q- Jacobi polynomials satisfy the orthogonality 
relations 



S ^Pn{q k )p m {q k ) w k (a,b) = 5 nm h n (a,b), 



k=0 



,k (h',Q)k (aq;q)oo u n 1-abq (q,bq;q) r 



(3.5) 



w k (a,b) = (aq) - — ' ' , h n (a,b) = - ' ' — r (ag) n 

(q; q) k (abq 1 ; q)^ 1 - abq n + 2 (aq, abq; q) n 

normalising ho(a,b) = 1. The little q- Jacobi polynomials satisfy 

L^p n (-;a,b;q) = X nPn (-,a,b;q), A„ = X„(a,b) = ,-(1 - 9 ")(1 - abq" +1 ) 

= (/w _ /W) + izl (/( £) _ /w) (3 ' 6) 

v / x x q 

In the context of Section 13.11 we take L 2 (p), respectively L 2 (i>), to be the weighted 
L 2 -space corresponding to the case (aq,bq), respectively (a, b). Note that 

w k (aq,bq)=r(q k )w k (a,b), r(x) = K~ l x(l-bqx), K = J^^H Z ^sj > °' ( 3 ' 7 ) 



In the context of Section 13.11 we see that we can give a five-term recursion formula for 
the operator T P,T defined by 

(T^f)(x) = ^(l- b qx)(bq 2 x-l){f(qx)-f(x)) 

1 p ( 3 - 8 ) 

+ — (1 - bqx)(x - l)(f(x/q) - f(x)) + a:(— (1 - 6gx) + r)/(x). 

In order to apply the link to 2 x 2-matrix-valued orthogonal polynomials we need to give 
the spectral decomposition of T P)T on L 2 (v) in another way. 

Proposition 3.4. Assume q^ 1 > b > 0. The operator T P,T for p = (1 + qyfab)(l + q 2 \fab), 
t = ^ (q\^ab(3+2q+bq 2 ) — bq(l+aq)} is a bounded self-adjoint operator with explicit spectral 
decomposition given by U: L 2 (v) — > L 2 (a) and UT = MU, where M is multiplication by 
2x/(^Jaq + (1 — K) / Jaq) and a is the normalized orthogonality measure for the continuous 
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dual q-Hahn polynomials [H] with parameters (A,B,C) = (^/qb,^/qb,q^/qb), and U is 
given by 

u ■■ . 9 = 7 x i-)> P n {-\ Vqb, Vqb, qVqb \ q) 

\Jw k {a,b) 

using the orthonormal polynomials on the right hand side. 

Remark 3.5. Recall that in [TT] we can introduce an additional degree of freedom, which 
is not possible in Proposition 13.41 On the other hand, considering more generally second- 
order difference operators on L 2 [y) we can introduce an additional degree of freedom in 
the parameters of the continuous dual Hahn polynomials, but then we have no longer a 
nice explicit expression for the 5-term recurrence as in Section [2j 



Proof. Let V: £ 2 (N) — > L 2 {y), e k i-> S q k / ^Jw k (a, b), be the unitary operator identifying 
the Hilbert space £ 2 (N) with standard orthonormal basis {ek}ken with the weighted L 2 {y) 
space for the little g-Jacobi polynomials. Let J P,T = V*T P,T V, then 

(y/aq + — ]= + -^L\j p,T e k = a k e k+1 + (3 k e k + a fc _ie fc _i 



Rxq \faq~' 
a k = (l-bq k+1 W(l-q*+i)(l-bq*+i), 



h = q k (bqV^q(l + q) + -±={1 + bq) - -£= + + q 2k (-^=(p - 1) - b 2 q 3 ^q) 



Comparing this Jacobi operator to the three-term recurrence relation for the orthonormal 
continuous dual g-Hahn polynomials with parameters (A, B, C) as in [TU §3.3], we see that 
we need {AB, AC, BC} = {bq, bq 2 , bq 2 } to get the right expression for a k . Since b ^ we 
get A = B, C = qB, and because of symmetry we obtain (A, B, C) = (\fbq~, \fbq~, qyfbq). 

In order to match the value of j3 k to the orthonormal continuous g-Hahn polynomials 
with these parameters we require 

bqy/Eq{\ + q) + -^(1 + bq)--^ + ^% = A + B + C + ABC, 
y/aq y?aq yfaq 



(p - 1) - b 2 q 3 ^ = ABC(1 + q- 1 ] 



/aq 

which determines the choice for p and r. Then (^Jaq + + J P,T has continuous 
spectrum [—2,2], which gives the statement on U. □ 

Now that we have determined for which values of (p, r) we have an explicit spectral 
decomposition in Proposition 13. 4} we have to work out the coefficients in Corollary 13.21 in 
this case. We start with ( 13. ip in this case, or equivalently 

p n (x; a, b; q) = a n , n p n (x; aq, bq; q) + a n ^ x Vn~\{x; aq, bq; q) + a nin _ 2 p„- 2 0; aq, bq; q). (3.9) 

By comparing leading coefficients in ( 13. 9 p we obtain 

(1 - ag n+1 )(l - abq n+1 ){\ - abq n+2 ) 



(1 - aq)(l - abq 2n+1 )(l - abq 2n+2 ) 



(3.10) 
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Using the orthogonality and (13. 7p we obtain 

oo 

a n -2,nh n -2{aq,bq) = } j p n (q k ; a, b; q)p n -2{q k ; aq, bq] q)r(q k )w k (a, b) 

k=0 

l n {a,b) 

using the expansion of p n -2{'] aq, bq; q)r(-) in terms of little g-Jacobi polynomials with 
parameters (a, b). This gives 

-bq n+2 (l-q n - 1 )(l-q n )(l-bq)(l-bq n ) 

a n ,n-2 - K ^ _ _ a69 3)(! _ afe^n-l)^ _ fl&g 2n) " ( " J 

Note that (13.111) is not clear from the general connection coefficient formula for little q- 
Jacobi polynomials due to Andrews and Askey, see [5j Ex. 1.33]. The coefficient a„ in _i 
can be obtained by comparing coefficients of x n ~ l on both sides. This gives, after a 
straightforward calculation, 

1 _ n (l-g")(l-q&^+ 1 ) ( 1 - aq n 1 - ag n+l \ 
° n ' n " 1 q (l-q)(l-aq) \l - abq 2 - 1 - abq 2n + 2 J 1 ' 

Using the orthonormal version we find that in this example the coefficients in (13. ip are 

q^ n (l - abq n+2 )^(l - abq 2 )(l - abq 3 )(l - abq n+1 ){\ - ag n+1 )(l - bq n+1 ) 
an ~ (1 - abq 2n + l ){\ - abq 2n + 2 )s/(l - abq n +±)(l- aq){l - bq) 

_i n _i - abq 2 )(l - abq 3 )(l - q n )(l - abq n+l ){\ - abq n+2 ) 

Pn — (7 2 G 2 , : 

(1 - q)y/(l - abq n+3 ){l - aq)(l - bq) 
1 — aq n 1 — aq n+l 



x 



1 — abq 2n 1 — abq 



2n+2 



(3.13) 



-bq 2n J (I - q n ~ l ){\ - q n ){l - aq){l - aq n ){l - bq)(l - bq r 

In ~ 



2ii^ 



aK (1 - a6g 2n - 1 )(l - abq 

Finally, we need the three-term recurrence relation for the orthonormal little g-Jacobi 
polynomials, which corresponds to (13. 2p with explicit values 

n ^/aq(l - ag n+1 )(l - &g" +1 )(l - g n+1 )(l - a6g"+ 1 )(l - abq n+2 ) 

On = q , 

(1 - abq 2n + l ){l - abq 2n + 2 )y/l-abq n + 3 , g u) 

t _ q n (1 - ag n+1 )(l - abq n+1 ) aq n (1 - q n )(l - bq n ) 



[1 - abq 2 ^ 1 )^ - abq 2n + 2 ) (1 - abq 2n )(l - abq 2n+1 ) 



We next want to use Theorem l2. ll with the spectral decomposition U given by Proposition 
13.41 so that we assume the situation of Proposition 13.41 The spectrum is simple, so that 
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VI2 = 0. It remains to calculate U<fto and U(j>\. Keeping track of normalization we have 



(TT , w ,s y/a{aq,bq,q; ffC,, 2 n 

[U(j) n )(cost) = — 7=== (og ;g) 



\/(a6g 2 ;9)c3o 

X ? "l 2 ' — Pk{cost] y/qb, y/qb, qy/qb I g) 

^ {q,bq 2 ;q) t 



k=0 



where we have used the standard notation, see [T4], for the continuous dual g-Hahn 
polynomials. Using one of the standard generating functions, see [TU (3.3.15)], and 
Pi(<l k ; a, b; q) = 1 - q k ^^ we find 

F (cost) = (£/0 o )(cost) = 



v/a(ag, fog, g; g)^ (fog 2 ; q)oo(qVb; q)^ ( y/qbe** , qy/qbe 1 * 

*M ..2 ^'V^e 



\/(abq 2 ;q)oo (fPy/QlQ)™ \ bq 

t? ( +\ (tta\< + \ V a ( a( l'bq,q;q) 00,, 2 , 

Pi (cost) = (f/0i)(cost) = —j===^(bq ;q) 00 ^ ^ 



x 



/ (qVb-^)^ ( y/qbe u , qy/qbe u _ 



1-afog) (g 2 A/fo;g)oo / y/qbe u , qy/qbe u s _ tt 

2^1 ( ,.,2 ^'^ 2e 



(1 - ag) (e^g^/gjg)^ \ & <? 

We summarize this situation in the following Proposition 13.61 using the explicit expression 
for the orthogonality measure cfer of the continuous dual g-Hahn polynomials, see [T4"j §3.3]. 



Proposition 3.6. Define the coefficients a n , b n and c n as in Corollaru \3.S\ with the explicit 
values for a n , (3 n , 7„ as in A3. 13|) . A n as (13. 6p . n , as in (13.14p and p ana 1 r as m 
Proposition \ 3.4\ The 2 x 2-matrix-valued orthogonal polynomials generated by the three- 
term recurrence relation (12.31) with initial conditions P_i(A) = 0, Po(A) = I satisfy the 
orthogonality relations 

(g, gfo, g 2 fo, g 2 fo; g)oo P (e ±2ii ;g)oo , r 
^ P n (cost)^ l( cost)P m (cost) (g6e±it>g26e±iSg26e±it;g)oo dt = 5 nm I 

with 

w (c ^ f) _ ( l^o(cost)| 2 F (cos*)Fi(cos*)\ 
KKllCOS J ~ \F (cost) F 1 (cost) |Pi(cost)| 2 )• 

In view of [TTJ §4] we view the 2 x 2-matrix-valued polynomials of Proposition 13.61 as 
matrix-valued analogue of (a subfamily) Askey- Wilson polynomials. The case b < can 
be dealt with similarly, where the case b = allows for additional degrees of freedom. 
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4. Example: spectral decomposition of an operator arising from quantum 

GROUPS 



In an influential paper [18] Koornwinder has introduced a special element p TjCr in the 
quantum £77(2) group. In this context it is important the have the action of this element 
in an infinite dimensional representation as an explicit 5-term recurrence relation. On the 
other hand, the spectral decomposition of the corresponding operator has been solved in 
[TT] exploiting the special case a — > oo as an intermediate step. So the spectral decompo- 
sition of the 5-term recurrence is completely known, and by the set-up of Theorem 12.11 we 
obtain orthogonality relations for 2 x 2-matrix-valued orthogonal polynomials with explicit 
coefficients for the three-term recurrence relation. The resulting Proposition 14.11 describes 
the weight function explicitly in terms of 2V?i-series. 

Throughout this section we assume a, r e R. In this case the Hilbert space is H = £ 2 (N) 
with standard orthonormal basis {f n }'^Lo- The operator T corresponds to the operator 
7r^(p T)Cr ) of [17, §6]; explicitly in the notation of ( 12.1 j) we have in this case 

a n = 1^/(1 - g 2n+2)( 1 _ ? 2n+4) > 

b n = -iq n+1 y/l - q 2n + 2 (e^(q- a - q a ) + e~^(g~ r - q T )) (4.1) 



Cm 



g 1+2n (cos(20) - l -{q~° - ff){q^ - <f)) 



Note that a_i = a_2 = fr-i = 0. Moreover, we have a symmetry (a, r, 0) (r, er, — 0) and 
(a, t, 0) -H- (— a, — r, + 7r). So we can assume a > r and a > —r. From [TTJ §6] we deduce 
that T has absolutely continuous spectrum [—1, 1] of multiplicity 2 and discrete spectrum 
(possibly empty) of multiplicity 1 at £_ US + , where, using the notation //(x) = i(x + x _1 ), 

£_ ={fi(-q 1 -' r -^ +2k ) fceN, g i— -+ 2fe > 1} 

(4.2) 

S + = { A t( ? 1 - CT+T+2fc ) | k e N, g 1 "^^ > 1}. v ; 

From [T71 §6] we can read off L 2 (V). Assume that a + r < 1, <r — r < 1, so that there 
is no discrete spectrum. Then V is a diagonal matrix with the orthonormal measure 
for the Al-Salam-Chihara polynomials with parameters (q 1+cr ~ T , — g 1_fJ_T ), respectively 
(q 1 ~ cr+T , — q 1+cr+T ), on the (1, l)-entry, respectively the (2, 2)-entry. Explicitly, /: [—1, 1] — > 
C 2 is in L 2 (V) if 

/»7T 

/ | /i (cost) | Vi (cost) + | / 2 (cos t)| 2 v 2 2 (cost) dt < oo (4.3) 
Jo 



r o, (g 2 ,-g 2 - 2r ;g 2 )oo( e ±2it ;g 2 )oo 



with 

wn(cost) = vn(cost; o r , o cr I a 2 ) , , , , 

V 7 U '^^ 2n(-q 2 -,q 2 ) O0 (q 1 +—e±* t ,-q 1 -—e±* t ;q 2 ) O0 (4.4) 

v 22 (cost) = Un(cost; q~ T , q~ a \ q 2 ) 
and so wi 2 (cost) = = t> 2 i(cost). 
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In order to write down the orthogonality measure for the 2 x 2-matrix-valued orthogonal 
polynomials from Theorem 12.11 we need to calculate Ufk for k = and k = 1 . Expanding 
the standard orthonormal basis into the basis {w^, ii^}^ =0 as in [T71 p. 410], and applying 
U, which is given by (Ai,A2) as in [TJ, p. 411], we get after a straightforward calculation 
that Uf k (X) = (UxMX), U 2 f k (X)y with 



m=0 iP-gSm 



OO 



. „ ,;-k p -ik4> ( 2r+2m\ 



(4.5) 



m=U II (jr 



where we have used the notation as in [TTJ Prop. 5.2, p. 410] for the length of the vector, 
the Al-Salam-Carlitz polynomials Pk{-) and the Al-Salam-Chihara polynomials h m (-). 
For k = we can use the generating function, see e.g. [T4"l (3.8.14)], directly to find 

F 1 , (cost;g r ,g CT | q 2 ) = (I7i/o)(cost) 

= 1 2fi 2 _ 2t ; g , ge 1 > 

(-g^g^ge'^g 2 )^ V -? / 

and ({7 2 /o)( cos £) = ^i,o ( c °st; g~ T , g _<J I q 2 ) is obtained from (C/i/ )(cos£) by replacing 
(a,r) by (-a, -r). 

For = 1 we have to take a linear combination. First, note, in the notation of [TTf 
Prop. 5.2],pi(x) = q~ T (l — g 2 )~ 5 (x + 1 — g 2r ), so that g 2m ) = — pi(q 2m+2,T ). In partic- 
ular, we obtain, also using [T7J Prop. 5.2], that (t/ 2 /i)(cosi:) is obtained from (U\f \){cost) 
by switching (a, r) to (—a, — r) and multiplying by —1. Using the same generating function 
for the Al-Salam-Chihara polynomials twice we obtain 

F lil (cost;g T ,g CT | q 2 ) = (U x f x ){cost) = 



>/(i - <f)(-9 2T ; ? 2 )oo H?vw);? 2 )oo ^ V -9 2 - 2T ; (4.7) 

2 i(2<f>-t) 



1 - q 2r ) {q 1+a ~ T e u , -q x - (J - T e it 



and (U 2 fi) (cos i) = — .^(cost; g~ T , g _<J | g 2 ). 

Proposition 4.1. Consider the matrix-valued polynomials P n generated by (12. 3p mt/i mi- 
tto/ conditions P_i(A) = 0, -Po(A) — -f where the entries of the matrices A n and B n 
are given by (14.11) wift a > r, o > —r. Assume moreover o + r < 1, a — r < 1 ; then the 
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matrix-valued polynomials P n satisfy the orthogonality relations 




P n (cost) W 2 (cos t)P m (cos t) * dt = S nm I, 



W 2 (cost)n = \F lfi ( y cost;q T ,q (J | q 2 )\ \ n {cost] q T , q a ) + ((a,r) ^ (-a,-r)) 

W 2 (cost) 21 = W 2 (cost) 12 = F lfi (cost;q T ,q a \ q 2 )v 11 {cost;q T ,q CT )F 1A (cost;q T ,q a | q 2 ) 

~ {far) (-a, -r)) 

W 2 (cost) 22 = \F 1A (cost;q T ,q a | q 2 ) | 2 u n (cos t; q T , q a ) + ((a,r) «-> (-a, -r)) 

where the functions on the right hand side are defined by (I4.6p . (14 .7p and the notation 
((a, r) <H- (—cr, — r)J means that we have to add the same term but with parameters (a, r) 
replaced by {—a, — r). 

Note that Videos is explicit as a sum of i + j terms, each term being a product of 
two 2V 9 i- ser i es - 

In case the assumption a + r < 1, a — r < 1 is dropped we obtain a finite discrete 
set of mass points in the orthogonality relations of Proposition I4.1[ and the weight W\ at 
these points can be calculated in the same way from Theorem 12.11 Alternatively, they can 
be obtained from writing the integral of Proposition 14.11 as a contour integral, and then 
shifting contours which leads to discrete masses at the poles with weights given in terms 
of residues analogous to the case of the Askey- Wilson polynomials, see pQ. 
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